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'Abstract 

\ Results of a long-term monitoring (> 10 years) of the broad line and continuum fluxes of three Active Galactic Nuclei 
C/3 (AGN), 3C 390.3, NGC 4151, and NGC 5548, are presented. We analyze the Ha and H/3 profile variations during the 
iy-j monitoring period and study different details (as bumps, absorption bands) which can indicate structural changes in 
the Broad Line Region (BLR). The BLR dimensions are estimated using the time lags between the continuum and the 
broad lines flux variations. We find that in the case of 3C 390.3 and NGC 5548 a disk geometry can explain both the 
■broad line profiles and their flux variations, while the BLR of NGC 4151 seems more complex and is probably composed 
of two or three kinematically different regions. 
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INTRODUCTION 



An important question in the study of Active Galac- 

j tic Nuclei (AGN) is the nature of the "central engine". 

Nowadays, it is widely accepted that the nuclear activity 
CN is caused by accretion of gas onto a supermassive black hole 
^ | (SMB H) . More than 40 years ago, Weedman fc Khachikvan 
'( 19681 ) classified AGN into two types: AGN1 whose spectra 
^ <3how broad (FWHM> 1000 km/s) and narrow (FWHM~ 
£N) 400 km/s) emission lines; and AGN2 with only narrow 
qq .emission lines. Later, the classification was improved by 
incl uding interm ediate types from AGN 1.2 to AGN 1.9 
Q\ '.(see lOsterbrockl Il989h . It was first believed that AGN1 
and AGN2 have different physical natures. However, more 
than 20 years ago a unified scheme of the AGN struc- 
ture was proposed, where an accretion disk (AD) around 
a SMBH are assumed to be in the center of all AGN. 
Around the very central part, two emitting regions, the 
"Broad Line Region (BLR)" and the " Nar row Line Region 
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angular size of the BLR is -Rblr < 0.001 arcsec, conse- 
quently even for nearby AGN1 this region is presently not 
resolved with the most sensitive optical interferometers. 
The only ways to study the BLR are indirect methods, as 
e.g. "reverberation mapping" (Peterson et al. 1993). 

It is well known that AGN vary in luminosity on time 
scales from years to hours, over the whole wavelength 
range from the radio to X-rays or 7-rays. In particu- 
lar, the flux in BELs varies with respect to the ionizing 
continuum flux with short time delays (days to weeks for 
Seyfert galaxies), due to light-travel time effects. If the 
BLR presents systematic motions such as infalling, out- 
flowing, circular motions, the profiles of the broad emis- 
sion lines should vary in a way related with its geometry 
and kinematics. The motions also depend on the processes 
of gas relaxation that follows changes in the ionizing flux. 

An important progress in understanding the BLR struc- 
ture was achieved as a result of multiwavelength monitor- 
ing campaigns performed within the framework of the In- 



(NLR)",are surrounded by a torus of dust (|Antonucci fc Mlld ternational AGN Watch a consortium organi zed to study 



19851 : Antonuccl . 1993f ). Then the difference in spectra be- 



tween AGN1 and AGN2 is related to the viewing angle of 
an observer with respect to the torus orientation. 

The BLR, i.e. the Broad Emission Lines (BELs) form- 
ing region, consists of gas obviously linked with the accre- 
tion process onto a SMBH. It is important to investigate 
its structure in order to use the BELs for the estimation 
of mass and the accretion rate of the SMBH. The typical 
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several Seyfert galaxies (IPeterson et al 

Its main results (see IPetersonl l2008l and references 
therein), and references therein) can be summarized as 
follows (i) The response time of the H/3 line to continuum 
variations varies from year to year and is correlated with 
the average continuum flux; (ii) the BLR sizes range from 
light-days to several light-weeks; (iii) the high ionization 
lines have a shorter response time than the low ionization 
lines, indicating the presence of an ionization stratification 
along the BLR; (iv) rotation or chaotic motions dominate 
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onto radial motions in the BLR; (v) the optical and ultra- 
violet continua vary almost simultaneously (i.e. without 
showing a significant delay). 

All these conclusions improved our knowledge of the 
BLR, but many questions remain still open, as e.g. the ex- 
istence of a small fraction of double peaked profiles which 
could indicate the presence an accretion disk, or that of 
an inflow or an outflow. To answer these questions, a long 
term monitoring of AGN with a detailed study of the dif- 
ferent broad line profiles was required. 

In this paper we give an overview of our recent investi- 
gation based on the long-time monitoring of three AGNs 
with different BEL profiles: 3C 390.3 with double peaked 
profiles, NGC 4151 with strongly variable BEL profiles, 
and NGC 5548 with shoulders in the profiles. 

2. OBSERVATIONS 

High quality spectra (S/N > 50) in the continuum 
around Ha and H/3 were obtained in the spectral range 
4000 to 7500 A, using a long-slit spectrograph equipped 
with CCDs, with a resolution between 5 and 15 A, us- 
ing the 6-m and 1-m SAO's telescopes (Russia) and the 
GHAO's 2. 1-m telescope (Cananea, Mexico). Spectropho- 
tometric standard stars were observed every nights. The 
image reduction process included bias subtraction, flat- 
field corrections, cosmic ray removal, 2D wavelength lin- 
earization, sky spectrum subtraction, addition of the spec- 
tra for each night, and relative flux calibration based on 
standard star observations. 

For the absolute calibration, the fluxes of the narrow 
emission lines are adopted for scaling the AGN spectra, 
because they are known to remain constant on time scales 
of tens of years (Peterson 1993). We thus assume that the 
flux of the [OIII] A 5007 line or [OI] A 6300 was constant dur- 
ing the monitoring period. Th e scaling of spectra was per- 
forme d by using th e method of[Van Groningen &: Wanders! 
( 19921 ) modified bv lShapovalova et alj (|2004N 1 I. This method 
allowed us to obtain a homogeneous set of spectra with 
the same wavelength calibration and the same [OIII]A5007 
flux. 

In order to investigate the long term spectral variabil- 
ity of AGN it is necessary to have a consistent set of spec- 
tra. Since Seyfert galaxies were observed with different 
telescopes, in different position angles, and with different 
apertures, first we had to perform corrections for the po- 
sition angle (PA), seeing and aperture effects. A detailed 
discussion on the nec essity for these corrections is given in 
Peterson et al.1 (1993), and will not be repeated here. 

Fluxes of 3C390. 3 and NGC 5548 w e re corrected only 
for aperture effects (jShapovalova et all l200ll 120041 ) . be- 



cause the host galaxy is weak and the nucleus have its 
angle size < 2" and spectra were taken with apertures 



4.2" x 19.8", 2.5" x 6", and fluxes were practically inde- 
pendent from seeing and position angle. The fluxes NGC 
4151 were corrected for the position angle (PA), seeing and 
aperture effects, because the host galaxy is very bright and 
brightness of nucleus is v ery big and its angle size > 2" (see 
Shapovalova et all 120081 for details). 

The mean error (uncertainty) for H/3 and the contin- 
uum around this line is < 3%, while it is ~ 5% for Ha. 
In order to study the broad components showing the main 
BLR characteristics, the narrow components of these lines 
and the forbidden lines were removed from the spectra. 
In this aim, we constructed spectral templates using the 
blue and red spectra in the minimum activity state. The 
broad and narrow components of H7, He II A4686, H/3 and 
Ha, were fitted with Gaussians. Then, we scaled the blue 
and red spe ctra according to our sca ling scheme (see the 



Appendix in lShapovalova et all 120011 ). using the template 



spectrum as a reference. The template spectrum and any 
observed spectrum were thus matched in wavelength, re- 
duced to the same resolution, and the template spectrum 
was subtracted from the observed one. More details about 



data reduction can be found in iShapovalova et al 
20081) . 
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1 see Appendix A in [Shapovalov a et al ] [|2004h 



3. RESULTS 

3.1. BEL profile variations in NGC 5548 

This nucleus was observed from 1996 to 2002. A de- 
tailed de scription of the observatio ns and results can be 
found in Sha povalova et al. I (120041) . The results can be 
outlined as follows: 

- The flux in the broad lines and in the continuum grad- 
ually decreased, reaching minimum values during May- 
June 2002. In the minimum state, the wings of Ha and 
H/3 became extremely weak, corresponding to a Syl.8 type 
and not to a Syl, as observed previously when the nucleus 
was brighter (see Fig.l). 

- The Ha and H/3 broad line profiles displayed double- 
peaked profiles, at radial velocities ~ ±1000 km/s with 
respect to the narrow components. The relative intensities 
and radial velocities of the peaks varied. During 1996, the 
red peak was the brightest, while in 1998 - 2002, the blue 
peak became brighter. Their radial velocities vary between 
- 500 and 1200 km/s (see Fig. 2). 

- In 2000 - 2002, a distinct peak was observed in the red 
wings of Ha and H/3, at a radial velocity of ~ +2500 km/s. 
Its radial velocity decreased from 2000 to 2002 by 500 — 
600 km/s (see Fig. 3). 

- On small time scales, the flux of different parts of 
the profiles changed in almost the same manner, being 
highly correlated both with each others (correlation co- 
efficient r ~ 0.94 — 0.98) and with the continuum flux 
(r ~ 0.88 — 0.97). These results indicate that the flux 
variability on short time scales is caused mainly by the 
reverberation effect (see Fig. 4). In Fig. 4 flux 1,2,3,4 
and shape 1,2,3,4 correspond to the flux or shape in the 
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~ 1998 Jun26; JD2450991 (high state) 
2002 Jun4; JD2452430 (low state) 
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Figure 1: The spectra of NGC 5548 corresponding to the high ac- 
tivity state (top) and the low activity state (bottom). Observed 
wavelengths (we recall that z=0.0167) are displayed on the X-axis, 
and fluxes (in units of 10 -15 erg cm -2 s — 1 A -1 ) are displayed on the 
Y-axis. 
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Figure 2: Year averaged observed profiles of the H/3 (respectively 
the Ho) broad emission component in NGC 5548, after subtraction 
of the continuum and the narrow components [OIIIJAA 4959, 5007 
(respectively [Nil] A6548, 6584 and [SII]A6717, 6731). The vertical 
lines correspond to the radial velocities ±1000) km/s - dashed lines, 
~ +2600 km/s - dash-doted lines, Okm/s - solid lines. 



3 



March-May 2000 




-10000 -5000 5000 

Velocity [km s~ ] 



Figure 3: Year-averaged normalized (with respect to the narrow H/3 
and [OIII] lines) profiles of H/3 of NGC 5548 in the period 2000- 
2002. The vertical dashed lines correspond to the radial velocity 
~ +2600 km/s. 



Ha and H/3 profile segments in the following radial velocity 
intervals: fluxl and shapel for (—3000, —2000) km/s; flux2 
and shape2 for (—1500, —500) km/s; flux3 and shape3 for 
(+500, +1500) km/s and flux4 and shape4 for (+2000, 



+3000 ) km/s. The shape determination is given in Shapovalova et al 
(2004, 3.3.1) using method suggested bv lWanders fc PetersonI 
(Il996h . 

- On long time scales, the variations of different parts 
of the line profiles were mildly correlated with each others 
and with the continuum variations, or simply did not cor- 
relate at all. So, the profile changes on long ti me scales are 

not du e to reverberation, as earlier discussed bv lWanders fc Peterson 
(Il996h . 

- Changes of the integral Balmer decrement were, on 
average, anticorrelated with the continuum flux variations. 
This is probably due to an increasing role of collisional 
excitation as the ionizing flux decreases. The behavior of 
the Balmer decrement in various parts of the line profiles 
was different in 1996-2000 as compared with 2001. 

- From CCF-analysis of the time delay between the 
continuum variations and the H/3 response during 13 years 
( Peterson et all l2002h . one finds that the H/3 time lag var- 
ied from of 6 days in 2000 to 26 days in 1998-1999. The 
relation between the size of the BLR and the ionizing lu- 
minosity is: i?(BLR) - L 05 . 

The variation in the broad line profiles of NGC 5548 
indicates presence of an accretion disk emission which con- 
tributes to the BEL fluxes. 

3.2. BEL profile variations in NGC 4151 

The galaxy was observed in period 1996-2006. The de 



tails about observations can be found in Shapovalova et al 



(|2008h . Here we outline some of the most important re- 
sults: 

- The continuum and line fluxes varied strongly (up to 
a factor 6) during the monitoring period. The emission 
was maximum in 1996-1998, and there were two minima, 
in 2001 and in 2005. As a consequence, the spectral type of 
the nucleus changed from a Syl.5 in the maximum activity 
state to a Syl.8 in the minimum state (Fig. 5). 

- There are three characteristic periods (1996-1999, 
2000-2001 and 2002-2006) in the profile variability. In 
the first period, when the lines were the most intense, 
a highly variable blue component was observed, which 
showed two peaks or shoulders at V r ~ —4000 km/s and 
V r ~ —2000 km/s in the rms Ha profiles and, to a less 
extent, in H/3. In the second period, the broad lines were 
much fainter; the feature at V r ~ —4000 km/s disappeared 
from the blue part of the rms profiles of both lines; only 
the shoulder at V r ~ —2000 km/s was present. A faint 
shoulder at V r ~ 3500 km/s was present in the red part 
of rms line profiles. In the third period a red feature (a 
bump or a shoulder) at V r ~ 2500 km/s was clearly seen 
in the red part of both the mean and the rms line profiles. 
The behavior of the rms profiles in the three periods indi- 
cates that the BLR has a complex structure and that its 
geometry probably changes with time. 
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Figure 4: Flux-flux, shape-shape and continuum flux-shape correlations for different segments of Ha and H/3 in NGC 5548 (the number of 
segments is noted inside the plot). The correlation coefficient r are given inside the plot. 
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Figure 7: The flux of different segments of the broad lines as a function of the continuum flux in NGC 4151. The line flux is given in 
erg cm _2 s -1 and the continuum flux in erg cm _2 s _1 A _1 
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Figure 8: The profile variations of Ho and H/3 in 3C390.3: observed (left side panels) and normalized (right side panels). Ha and H/3 are 
denoted with solid and dashed line, respectively. 
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- 2005 Jun.8 JD2453530 (low state, GHAO's 2.1m telescope) 
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Figure 5: The spectra of NGC 4151 corresponding to the high ac- 
tivity state (top) and to the low activity state (bottom). Observed 
wavelengths (we recall that z=0.0033) are displayed on the X-axis, 
and fluxes (in units of 10 -14 erg cm -2 s _1 A -1 ) are displayed on the 
Y-axis. 
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Figure 6: Some examples of month-averaged profiles of the Ha and 
H/3 broad emission lines from 2000 to 2006 in NGC 4151. The ab- 
scissae shows the radial velocities relative to the narrow components 
of Ha or H/3. The vertical lines correspond to radial velocities: 
— 2600 km/s; Okm/s and 3000 km/s. The profiles are shifted ver- 
tically by a constant value. 



In 2005 (May- June), when the nucleus of NGC 4151 
was in the minimum activity state, the line profiles had a 
structure with two distinct peaks (bumps) at radial veloci- 
ties of (-2600, +2000) km/s in H/3 and (-1300, +2300) km/s 
in Ha. In 2002 a distinct bump appeared in the red wing 
of both lines at +3100 km/s. The radial velocity of the 
bump in the red wing changed from +3100 km/s in 2002 
to- 2100 km/s in 2006 (Fig. 6). In 1999-2001 we observed 
a broad deep absorption line in H/3 at the radial velocity 
V r ~ -400 km/s. 

- The CCFs are often asymmetrical and the time lags 
between the lines and the continuum are badly defined, 
indicating the presence of a complex BLR, with dimensions 
of 1 to 50 light-days. 

- The Balmer decrement (BD = j^) was maximum in 
1999 — 2001. An anti-correlation with the continuum was 
observed in the form of two almost parallel scries corre- 
sponding to 1999 - 2001 and 2002 - 2006 with a difference 
in BD by a factor - 1.5. In 1996 - 1998 the BD did 
not show any dependence with the continuum variations. 
The fact that the BD had the different values during the 
monitoring period (as well as different values along the 
profiles) also indicates a multi- component origin of the 
broad lines. It may be caused either by absorption, or by 
different physical condition s in different parts of the BLR 
(jShapovalova et al.l . l2009al ). 

- There was a linear relationship between the emission 
line and the continuum flux variations when the contin- 
uum flux Fc was small (Fc < 7 • 10~ 14 erg cm _2 s _1 A _1 ). 
The fact that the Ha and H/3 fluxes were well correlated 
to the continuum flux means that the ionizing continuum 
was a good extrapolation of the optical continuum. When 
Fc > 7- 10~ 14 erg cm~ 2 s _1 A _1 , the line fluxes either were 
weakly correlated, or simply did not correlate at all with 
the continuum flux. So it seems that the lines were satu- 
rated at high fluxes (see Fig. 7). 

3.3. BEL profile variations in 3C 390.3 

The observations of 3C 390.3 were performed from 1995 
to 2007. The results of t he first part of the monit oring 
period were publi s hed i n IShapovalova et alJ (|200ll ) and 



Shapovalova et all d2008l) 



be published soon ( S hapovalova et al 



and those of the la st part will 
" 2009b[) . The main 



results can be outlined as follows: 

- The broad component of the H/3 and Ha lines and the 
continuum flux varied by a factor of 3-4 during 1995-2007. 

- The H/3 profiles changed drastically: the wings were 
very strong in 1996-1998 (like those of a Syl), with a 
prominent bump in blue wing, at Vr — — 3700km/s with 
respect to the narrow H/3 component, and with a less in- 
tense flat bump in the red wing; in 1997 the wings were 
very weak (like those of a Syl. 8). The blue wing was al- 
ways brighter than the red one in 1995-2007 (see Fig. 8). 

- The shift of the blue peak in the integral H/3 profile s 
closely follows the trend found by lEracleous et aT ( 19971 ). 
These authors have shown that, if the displacement of the 
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Figure 9: The Auto- and Cross-Correlation Functions in 3C390.3. 
In panel (a), the narrow curve represents the sampling window for 
the autocorrelation function (ACFW) between the spectral contin- 
uum and the H/3 line. Also, the autocorrelation functions (ACF) for 
H/3 (1) - and the continuum light curves (2) are shown. On panels 
(b) and (d), the CCFs for the periods 1995-1999 and 1997-1999 are 
shown. Curves (1) represent the CCF between the H/3 flux and the 
continuum at 5125 A; curves (2) represent the CCF for the H/3 flux 
and the continuum, and label (3) refers to the CCF computed as 
in case (2), but for a restricted interpolation of 100 days. Panel (c) 
shows the CCF between: (1) - H/3 blue and red wings light curves, 
(2) - H/3 red wing and H/3 core fluxes and (3) - H/3 blue wing and 
H/3 core light curves. 



blue bump peak velocity originates in individual regions 
associated with a massive binary BH, then the inferred 
rotation period would be about 800 years, and the corre- 
sponding mass would then be larger than 1O 11 M . They 
argued that such large binary black hole masses are diffi- 
cult to reconcile with other observations and with theory. 
Therefore, they reject the hypothesis of a binary black hole 
for 3C 390.3. Our results provide further support to their 
conclusion. 

- From the cross-correlation analysis, the time lag of 
the H/3 emission line response relative to the continuum 
variations was found to vary with time: during 1995 — 1997 
a time lag of about 100 days was detected, whereas for the 
period 1998 — 1999 a double-valued time lag of ~ 100 days 
and ~ 35 days (double-peaks in the CCF) is present in our 
data (see Fig. 9). 

The observed flux of the H/3 wings varied quasi- 
simultaneously with the same lag relative to the continuum 
variations. Also, the CCF analysis of the H/3 wings did not 
reveal any delay between the variations of the wings and 
the central part of the line, or between two parts of the 
wings. This behavior is a strong evidence for predomi- 
nantly circular motions in the BLR. 

- From the profile differences of H/3, we found that 
the variations in velocity of the blue and red bumps and 
their differences were anticorrelated with the flux varia- 



tions of the H/3 and with its time lag with respect to the 
continuum. It means that the maximum of the line emis- 
sivity moves across the disk and corresponds to smaller 
radii when the continuum flux decreases (the bump veloc- 
ities increase) , and to larger radii when the continuum flux 
increases (the bump velocities decrease). These transient 
phenomena are expected to res ult from a variable rate o f 
accretion close to the black hole. IShapovalova et al.l (|200ll) 
concluded that these results favor the formation of the H/3 
broad emission line in the accretion disk. The modeling of 
the H/3 integral profiles was performed within the frame- 
work of the accretion disk model, and a satisfactory agree- 
ment for the disk with an inclination of 25°, and with a 
region of maximum emission located at about 200i?g was 
obtained (R g : gravitation radius) . 

Recentlv lStalevski et all (j2009l ) fitted the profil e of H/3 



in 3C 390.3 (1995-1999) using a perturbation model (St alevski et al 
20081) and found that S-like perturbation in the disk emis- 
sivity can explain the observed line profile variations. 

4. CONCLUSIONS 

Here we briefly outline some results obtained using 
long-term monitoring campaigns of the three AGN. We 
investigate the variations of the H/3 and Ha line fluxes 
and profiles, as well as the correlations between the line 
and the continuum flux variations. Our conclusions for 
the three objects are: 

1. For NGC 5548 our results favor the formation of 
the broad Balmcr lines in a turbulent accretion disk con- 
taining big moving optically thick inhomogeneities, able to 
reprocess the central continuum. 

2. The BLR of NGC 4151 seems much more complex. 
The broad lines are likely produced at least in two or three 
distinct regions with different physics (one photoionized 
by the AGN central source and another non-radiatively 
heated). 

3. In 3C 390.3 the double-peaked broad emission lines 
clearly indicate a disk emission. The variations in line 
profiles can be explained as a perturbation in an accretion 
disk (more details will be given in forthcoming paper by 



S hap ova lova et al 



(|2009bl) ). 

In Table 1 we compare some characteristics of varia- 
tions of the three AGN. As can be seen, the variability 
can be strong (factor 3-6 in the spectral lines), and dur- 
ing a long period the type of an AGN can change, e.g. 
from a Syl to a Syl.8. Two objects show evidence of disk 
emission. 

It is interesting that during a long period, not only the 
dimension, but also the geometry of the BLR, can change 
(as in NGC4151). All these results should be taken into 
account when BELs are used to estimate the SMBH masses 
and accretion rates. 
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